Preprint typeset in JINST style - HYPER VERSION 



Isotopically modified Ge detectors for Gerda: from 
production to operation 



D. Budjas , M. Agostini', L. Baudis^ E. Bellotti'^, L. Bezrukov^, R. Brugnera^, C. 
Cattadori^, A. di Vacri", R. Falkenstein^, A. Garfagnini^, S. Georgi ', P. Grabmayr*^, A. 
Hegai^ S. Hemmer^ M. Huir, J. Janicsko Csathy', V. Kornoukhov^'', B. Lehnert*, A. 
Lubashevskiy'', S. Nisi G. Pivato^, S. Sclionert', M. Tarka', and K. von Sturm*^ 

" INFN Labomtori Nazionali del Gran Sasso LNGS, Assergi, Italy 

^ Institut fur Kern- und Teilchenphysik Technische Universitdt Dresden, Dresden, Germany 

Institute for Reference Materials and Measurements, Geel, Belgium 

Max Planck Institut fiir Kernphysik, Heidelberg, Germany 

Universita di Milano Bicocca, Milano, Italy 
^ INFN Milano Bicocca, Milano, Italy 

* Institute for Nuclear Research of the Russian Academy of Sciences, Moscow, Russia 

Institute for Theoretical and Experimental Physics, Moscow, Russia 
' Physik Department El 5 and Excellence Cluster Universe, T. U. Milnchen, Germany 
^ Dipartimento di Fisica e Astronomia dell' Universita di Padova e INFN, Padova, Italy 

Physikalisches Institut, Eberhard Karls Universitdt Tubingen, Tubingen, Germany 
' Physik Institut der Universitdt ZUrich, ZUrich, Switzerland 
E-mail: |qrabmayr@uni-tuebinqen . de| 

Abstract: The Gerda experiment searches for the neutrinoless double beta (OvjSjS) decay of 
^^Ge using high-purity germanium detectors made of material enriched in ^^Ge. For Phase II of 
the experiment a sensitivity for the half life T^^^ ~ 2 • 10^^ yr is envisioned. Modified Broad En- 
ergy Germanium detectors (BEGe) with thick n+ electrodes provide the capability to efficiently 
identify and reject background events, while keeping a large acceptance for the OvjSjS -decay signal 
through novel pulse-shape discrimination (PSD) techniques. The viability of producing thick- 
window BEGe-type detectors for the Gerda experiment is demonstrated by testing all the produc- 
tion steps from the procurement of isotopically modified germanium up to working BEGe detectors. 
Comprehensive testing of the spectroscopic as well as PSD performance of the Gerda Phase II 
prototype BEGe detectors proved that the properties of these detectors ai^e identical to those pro- 
duced previously from natural germanium material following the standard production line of the 
manufacturer. 

Furthermore, the production of BEGe detectors from a limited amount of isotopically modified 
germanium served to optimize the production, in order to maximize the overall detector mass yield. 
The results of this test campaign provided direct input for the subsequent production of the enriched 
germanium detectors. 
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1. Introduction 

Gerda []I], 0] is an experiment located at Laboratori Nazionah del Gran Sasso (LNGS) of INFN 
that searches for the neutrinoless double beta (Ovj8j3) decay employing high-purity germanium 
(HPGe) detectors enriched in the isotope ^^Ge. The bare detectors are operated in liquid argon in 
order to minimize the background from radioactive isotopes in materials near the detectors. The 
first phase of the experiment is employing reprocessed detectors (18 kg total mass, enriched to 
86% in ^""Ge) from the past experiments HdM and IGEX [0, §]. Gerda aims at a reduction of 
the background in the region of interest at Qpji= 2039 keV to a level of 10^^ cts/(keV-kg-yr), 
which is about a factor of 10 lower than previous state-of-art. So far, a background index (BI) of 
~2-10^^ cts/(keV-kg-yr) has been reached [Q]. In the absence of a signal and given the current 
BI, Gerda expects to set 90 % probability lower limits of Ti/2 > 1.9-10^^ yr for an exposure of 
20 kg-yr in Phase I [0] to scrutinize the claim In order to improve the half-life sensitivity by a 
factor 10 in Phase II of the experiment, 20 kg of new detectors from enriched germanium ('^"'Ge) 
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will be added to reach the proposed exposure of 100 kg-yr within reasonable time. Simultaneously, 
the BI is expected to be reduced to 10^^ cts/(keV-kg-yr) or below lU]. 

Apart from improved shielding, active background suppression techniques will be employed 
to reach the envisioned BI of ~10^^ cts/(keV-kg-yr). The measures comprise the instrumentation 
of the liquid argon and new enriched germanium detectors. The chosen detector type for Phase II 
is a modified thick window Broad Energy Germanium (BEGe) detector manufactured by Canbena 
Olen, Belgium [Q]- Their proper functioning in LAr has been demonstrated [0]. In comparison with 
the presently used coaxial germanium detectors, BEGe detectors have lower noise which translates 
into better energy resolution. Furthermore, their characteristic internal electric field distribution al- 
lows the identification of OVj8j8-decay like signals and rejection of background signals with higher 
sensitivity and robustness using a novel technique of pulse-shape discrimination (PSD) [0, 0]. 

Compared to coaxial HPGe detectors used in Gerda Phase I and past Ovj3j8 experiments [EL 
m , the BEGe detector design imposes tighter constraints on the impurity and defect concentrations 
in the germanium crystal material. These requirements arise due to the non-coaxial electrode ar- 
rangement in BEGe detectors and the resulting electric field profile; these are more demanding to 
achieve complete charge collection within the detector volume. 

To demonstrate that working BEGe detectors with the desired properties can be produced from 
the procured '^"Ge material, a comprehensive validation was performed. The tests of the production 
chain started from procurement of isotopically modified germanium material of the same chemical 
history as the already available '^"Ge, i.e. the leftovers from the enrichment process, depleted 
in ^^Ge, that is named '^'^PGe. Further steps included the purification, crystal pulling and diode 
production. The tests finished with spectroscopic and PSD performance characterizations of the 
manufactured detectors. 

A second goal was the optimization of the production processes in order to maximize the mass 
yield of functional detectors from the available amount of '^"'^Ge material. These improvements 
were directly implemented into the production process of the Gerda Phase II '^"'^Ge detectors. 
The data obtained with the characterization of the '''^PGe BEGes provide an important reference for 
assessing the performance of new '^"'^Ge BEGes. 

The various steps of the production chain are described in the following section 0. Section 
reports the results of the detector characterization. 

2. Procurement and production 

The supply chain starts with the procurement of isotopically modified germanium in the form of 
Ge02. The isotopic composition is then confirmed by different types of material analyses per- 
formed at different laboratories. The production process continues by reduction of the oxide to its 
metallic form and subsequent purification via zone refinement. This is followed by crystal pulling, 
the determination of electrically active impurity concentrations in the crystal volume, and selection 
and cutting of slices useful for detector production. Detectors are fabricated by transforming the 
germanium crystal slices into junctions that are then mounted into vacuum cryostats for tests of 
their properties. 
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2.1 Material procurement and isotopic analysis 

In total, 34 kg of '^^PGeOi material was acquired in early 2009 from ECP Zelenogorsk, Russia mm . 
This material was a by-product from the previous enrichment process for Gerda and was chem- 
ically processed identically as the enriched germanium material. Given the higher abundance of 
'''^PGe after the enrichment process with respect to the '^"''Ge, the cost per mass of '''^PGe is sub- 
stantially lower. The only physical difference of the chemical identical material, is the isotopic 
composition and in particular the abundance of the isotope ^^Ge. 

Several 1 g sized samples of the Ge02 material were used for isotopic analysis via inductively- 
coupled plasma mass spectrometry (ICPMS) and neutron activation analysis (NAA). The combined 
results for the isotopic abundance fce, i-e- the relative number of nuclei, are shown in Table |T]. The 
standard deviation is dominated by the different results from the various methods rather than by the 
individual statistical uncertainties of the measurements. An isotopic analysis of the material after 
completing the full production chain is envisioned. Such a test could detect possible admixtures of 
natural germanium ("^'Ge) during the production processes. In order to maximize the sensitivity 
for such tests, the isotopic abundance of ^^Ge in the procured material was required to be <1 % 
from the beginning. 

Table 1. Isotopic abundance /ce of the depleted Ge02 samples, averaged over all measurements performed 
at LNGS (ICPMS), INR RAS, Moscow (ICPMS), University of Tubingen (NAA, with sample irradiation at 
FRM II Munich) and IRMM, Geel, Belgium (NAA, with sample irradiation at SCK»CEN, Mol, Belgium). 
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fce 


0.223(8) 


0.300(4) 


0.083(2) 


0.388(6) 


0.006(2) 



2.2 Crystal production 

After taking samples for ICPMS the input was 33.84 kg of GeOi material that was reduced to its 
metallic form and purified via zone refinement at PPM Pure Metals in Langelsheim, Germany WHW . 
This process yielded 21.49 kg of 6N grade germanium with specific resistivity Pvpec >50 Hem 
and ~1.37 kg beyond this specification. Using the specified abundances and the atomic mass 
Ma=72.44 the input mass amounts to 23.47 kg metallic germanium. Thus, the purification process 
was performed with an efficiency of 91.6 %. 

Crystal pulling and detector slice production was performed at Canberra Oak Ridge, USA 013]. 
The manufacturer modified its production process with the objective to maximize the total mass 
yield of final detectors for a fixed amount of starting germanium material. In total four p-type 
crystals were produced from which 14 detector slices were cut with a total mass of lO.O kg. This 
corresponds to a yield of 46.5 % with respect to the starting 21.49 kg of 6N depleted germanium 
material. Given the exploratory goal of this production test, also non-standard crystal slice dimen- 
sions and shapes were tested. Moreover, slices close to the seed-end and tail-end side of the crystal 
were investigated, to test whether working detectors can be produced. The number of test and 
inspection slices, usually I mm thick, was minimized. 
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Figure 1. Summary plot of the '^''PGe crystals and their slices. The orientation of the crystals is such that 
the seed end is on the top. The slices processed into detectors are shown in green, and the one slice which 
was turned into a functional diode, but not mounted into a vacuum cryostat is shown in blue. The thicker 
border line of the processed slices indicates the approximate extent of the Li-infused layer. The masses of 
the slices are given in g. 



Prior to defining tlie dimension of each detector slice, a full electrical field calculation was 
carried out for the conceived detector geometry, redundantly by the manufacturer and by us. For 
the latter the method of field calculations is summarized in Ref. O. Design criteria required that 
the depletion voltage should not exceed 4 kV and that all regions inside the detector volume should 
have sufficient strong electrical fields to ensure full charge collection. The diameter of the final 
slices vary from 72 mm to 78 mm, the height from 23 mm to 41 mm and the mass from 350 g to 
930 g. A schematic view of the crystal slices is depicted in Fig. [1]. The detector slice identification 
used in this report follows the order of crystal production and a letter corresponding to one of 
maximal four slices cut from each crystal; with A being the seed-end slice and D the tail-end slice. 

2.3 Detector manufacturing 

The detector manufacturing was finalized at Canberra Olen, Belgium [Q], by creating the Li-infused 
n+ electrode covering most of the surface, a B-implanted small-area p+ electrode located in the 
middle of the bottom flat face of the cylindrical slice and a circular non-conductive groove sep- 
arating the two electrodes. The modified BEGe detectors for Gerda differ from the standard 
production version [1T3I1 . The latter have a thin entrance window on the top flat surface to increase 
the sensitivity to low energy photons, while the GERDA version has a 0.4 mm to 0.8 mm Li-infused 
n+ electrode also on this part of the detector surface. 

The detectors were mounted in standard vacuum cryostats, equipped with Canberra 2002CSL 
charge sensitive preamplifiers. Five crystal slices were transformed into detectors with at least 
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one slice taken from each crystal. The mass loss turning a crystal slice into a diode is usually of 
the order of ten grams and thus negligible. However, slice 2/A was turned into a working diode 
only after removal of approximately 160 g of material from the seed-end side. This diode was not 
mounted in a vacuum cryostat because of its chamfered shape, which is incompatible with standard 
detector holders. 

The mass yield of crystal slices can be further increased to approximately 55 % employing 
twice the initial mass, as employed in the production of enriched detectors. Residual crystal pieces, 
that are not adequate for detector production, can return into the zone refinement or crystal pulling 
process. Moreover, it is estimated that a large fraction of kerf and grinding losses can be recovered, 
chemically reprocessed and returned to the zone refinement and crystal pulling process. 

2.4 The time line 

The procurement and production has been achieved in a rather short time. It proved rather helpful 
that: (/) negotiations on buying isotopically modified material form ECP started in Moscow on 
February 26, 2009. The material arrived in Tiibingen on April 30. Within one week 40 samples 
were analyzed regarding the isotopic composition by four institutes. (//) Zone refinement at PPM 
Pure Metals was accomplished by June 2009. (///) The 6N material arrived at Oak Ridge on August 
7. {iv) The first shipment of three detector grade slices left by plane for Olen on December 16, 2009. 
(v) The first two detectors were ready on March 19, 2010. (v/) The second batch of two detector 
slices was sent to Olen on June 15, 2010. {vii) The last shipment occurred on January 11, 201 1. 

The larger periods between the shipments was owing to intermediate testing at Laboratori 
Nazionali del Gran Sasso. Thus, within less than 2 years the production chain could be setup 
and verified. Having understood the details, the enriched BEGe detectors could be produced even 
faster. The main difference concerns a more carefully shielded transport to avoid activation of the 
enriched material. 

3. Detector characterization 

The five detectors produced underwent a comprehensive testing campaign to characterize their 
charge collection, spectroscopic and PSD performance, as well as long-term stability. An analog 
data acquisition (DAQ) chain consisting of a standard spectroscopy amplifier set to 10 /is shaping 
time and a multi-channel analyzer was used for the spectroscopic measurements. Digital recording 
of pulse shapes for subsequent off-line PSD analysis was accomplished using a DAQ system con- 
sisting of a non-shaping amplifier and a Struck SIS 3301 14-bit flash analog-digital converter with 
a read-out speed of 10^ samples per second similar to the one used for the Gerda acquisition at 
LNGS [0]. The recorded pulse-shape data were processed and analyzed using GELATIO finH . a 
Gerda software framework for advanced data analysis and digital signal processing. 

Four out of five detectors show a rather similar behavior resulting in the same parameters 
comparable to the standard detector. Only detector 1/D exhibits a reduced performance, the origin 
of which is not fully clarified. It was the last detector delivered and was tested separately. Various 
conjectures will be discussed in the following, e.g. noise contributions or inhomogeneous impurity 
distribution. However, none of these hypotheses can be finally rejected or proven. 
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3.1 Operational characteristics of the detectors 

The first set of measurements was carried out to investigate the operational characteristics of the 
detectors. This included the determination of: (i) the depletion voltage, (ii) the dead layer, (Hi) 
the active volume and other active dimensions, and (iv) the uniformity of charge collection along 
detector surfaces. 

The tests were iterated to ensure that a stable optimum for the parameters had been reached. 
Table summarizes the main results along with the equivalent results for a reference detector - 
a BEGe from standard Canberra production made of natural germanium, with 81 mm diameter, 
32 mm thickness and a mass of 878 g QBD- The particular's of the individual measurements ai^e 
discussed below. More details about the dead layer and active volume measurements can be found 
in chapter 8 of Ref. O. 

Particularly useful are scans along the surfaces of the detectors. A ^^^Am source was coUi- 
mated to provide a 59 keV photon beam within a spot with a diameter of 1 mm on the detec- 
tor surface. Due to the coUimation, only a small fraction of emitted photons reach the detector. 
For ease of comparison between different detectors and between results and Monte Carlo predic- 
tions the normalized counts rates /?„, namely detected photons over emitted photons, are plotted 
as function of position. The activity of the sources is known to 1 .5 % uncertainty. The intensity 
distributions /?„ will be compared with Monte Carlo (MC) simulations as shown in Fig. 0. 

3.1.1 The depletion voltage 

The reverse-bias voltage at which the detector volume is maximally depleted of free charge car- 
riers, referred to as the depletion voltage V^ep, was determined first. The corresponding volume 
in the detector is called 'active volume' Vact, while the remaining part, in general the surface with 
the 'dead layer' and electrodes, remains inactive. The depletion voltage is found by irradiating 
the detector with a ^°Co source while scanning the bias high-voltage range and determining the 
minimal voltage at which the spectroscopic characteristics (rate, width and position of 7 lines) are 
unchanged. Table ^ also includes the operational bias voltage Vop as recommended by the man- 
ufacturer. Vop is a conservative value above the full depletion clearing the active volume free of 
charge caiTiers. The excess of Vop over Vjep determined by the bias voltage scan is typically 1 kV. 

3.1.2 The dead layer 

To determine the thickness of the Li-infused n+ surface dead layer on the top area of the detector, 
the intensities of ^^^Am and ^^^Ba low-energy / lines are compared with Monte Carlo (MC) simu- 
lations as described in ref. [ITSI . |T3, 0, [TEIl. The values of the measured dead layer, dd, are collected 
in Table 0. The uncertainties include, besides the statistical and fit uncertainties, also a measure of 
the systematic uncertainty of the determination methods, expressed as the standard deviation of the 
results obtained with the ^^^Am and '^^Ba methods. Knowing the thickness dd the involved mass 
mrf can be calculated. A first estimate for the active mass fraction from this geometrical consider- 
ation can made: fg^o = {mtot —>nd)/mtot- Values in the range (91 - 96) % are obtained from this 
procedure as shown in Table 0. 
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Figure 2. Normalized intensities /?„ for the 59 keV peak of the ^"^^ Am source compared to MC predictions. 
The respective insets show the detector vacuum housings and holder cross sections with the beam position 
and direction indicated. Left: R^ y from the horizontal scans in x- and y- direction for the detector 4/D. Right: 
from a vertical scan along the side surface of detector 4/C. 



3.1.3 The active volume 

The active volume Vact is determined by irradiating the detector with ^"Co and reproducing the 
1332.5 keV y line detection efficiency with MC simulation while assuming that the only inactive 
part of the detector volume is the Li-infused n+ surface dead layer (c.f. ref. [1T51 IT^ . ITTH ). The 
uncertainty is dominated by the uncertainty of the source intensity and of its position with respect 
to the detector. In determining the active fraction of the detector volume /„,,, a small contribution to 
its uncertainty arises also from the uncertainties of the outer geometrical dimensions provided by 
the manufacturer. Comparing the results in Table 0, the nice agreement of fav with fgeo is observed. 
Only detector 1/D shows a strong deviation that can not be explained so far. 

3.1.4 Determination of uniformity from surface scans 

The characterization of the uniformity of charge collection along the detector surfaces is performed 
by a scan, that observes the variation of spectroscopic parameters (peak position and energy resolu- 
tion) along the scanning points. The variation of these parameters was found to be smaller than the 
random fluctuations caused by electronic disturbances in the DAQ setup: typically below ±0.1 keV 
for the 59.5 keV peak position and below ±50 eV for FWHM. 

Secondly, the uniformity of the dead layer thickness was measured by observing variation of 
the peak count rate along the scanning path. The uniformity along the top and side surfaces is con- 
firmed within the Umits of the measurements. This method can provide also an additional measure 
of the absolute thickness of the dead layer. However, the accuracy is much reduced compared to 
the methods described in section |3.1.2| due to the added uncertainties of the source activity and 
coUimation; the latter being dominant in the order of 10 %. 
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Table 2. Compilation of selected operational characteristics determined for the prototype detectors includ- 
ing the operational voltage Vop given by the manufacturer, the depletion voltage Vjep, dead layer thickness 
dd, active volume fraction calculated from the dead layer fgeo and fav as determined by the ^'^Co scan, 
the active diameter 0act, and the active height haci- For comparison the parameters of a reference natural 
germanium BEGe from the standard Canberra production are shown. 



detector 


Vop 


Vdep 


dd 


fgeo 


fav 


0act 


hact 


crystal/slice 


[kV] 


[kV] 


[mm] 


[m^/m^] 


[m^/m^] 


[mm] 


[rmn] 


reference 


4.5 


3.7 


0.430(23) 


0.9534(25) 


0.949(18) 


80.0(3) 


31.0(4) 


1/D 


3.5 


2.5 


0.348(15) 


0.9530(20) 


0.874(13) 






2/B 


4.5 


3.5 


0.761(15) 


0.9146(16) 


0.911(12) 


74.3(6) 


29.8(4) 


3/B 


5.0 


3.7 


0.43(4) 


0.957(4) 


0.956(15) 


73.5(5) 


39.5(4) 


4/C 


4.0 


3.0 


0.64(5) 


0.931(5) 


0.919(13) 


72.2(9) 


31.8(5) 


4/D 


4.5 


3.5 


0.394(12) 


0.9544(14) 


0.932(15) 


72.7(9) 


31.7(5) 



3.1.5 Determination of active lengths by surface scans 

The surface scans can be evaluated furthermore in order to extract the active diameter 0act at the 
very top and the active height hact- From both parameters an independent measure for the active 
volume Vact can be derived in principle. However, since some detectors have a conical shape as 
shown in Fig. |I], the evaluation of Vact of these detectors can not be accurate. The active lengths are 
measured by detecting the reduction of the peak count rate at the edges of the active volume in the 
scans. The obtained values reported in Table have relatively large uncertainties due to the limited 
positioning accuracy of the scanning device and the relatively large size of the scanning steps. 

The intensity distributions and Ry resulting from the scans on the top of the diode along the 
X- and j-directions, respectively, are compared with MC simulations as shown in the left part of 
Fig. ^. The rates are low due to collimation and thus multiplied by 10^ for display. The uniformity 
of the charge collection along the top is confirmed. As can be noticed, only some minor discrepan- 
cies between the simulations and the measurement results were observed, particularly at the edges. 
These are most likely to be attributed to an imperfect aUgnment of the scanning device, which was 
accurate to about 1 mm. 

The right part of Fig. @ displays R^ obtained from the scan along the side surface of detector 
4/C. The visible variations in count rates are due to blocking of the coUimated beam by the copper 
holder of the detector, which has a non-uniform shape. The top of the detector protrudes by ~1 mm 
from the copper holder, around z = 6 mm. Due to the ~0.5 mm dead layer, the exposed active part 
of the detector at this location is significantly smaller than the ~ 1 mm diameter of the coUimated 
beam. Therefore the count rate at that position is particularly sensitive to the vertical position of 
the detector in the holder. Small inaccuracies in the MC model thus lead to the visible discrepancy 
between the measurements and the MC data. 

3.1.6 Summary of active and dead volume determinations 

The results of these measurements are summarized in Table 0. The measurements of the geometry 
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of the active volume are in general consistent with the expectation that the full volume of the 
detector is active except for the surface dead layer. Furthermore, the determinations of the dead 
layer using the comparison of the rate in the collimated scans with MC simulations were also in 
agreement with the other methods, despite their lower accuracy. 

The detector 1/D has been an anomaly in this respect: the measured surface dead layer makes 
up (4.7=b0.2) % of the total detector volume while the active fraction fav measured with ^°Co 
indicates that only (87ibl) % of the total volume is active; compare also fgeo- The deviation could 
be an indication of a deficiency in charge collection inside the detector volume. Due to the limited 
accuracy, the scans along detector surfaces were inconclusive with regards to a possible dead layer 
thickness variation. Nevertheless, this discrepancy is not large enough to pose a serious concern 
for the detector production for the enriched BEGe detectors of Gerda. 

Actually, as a result of these tests the design of the detector holder was changed from the 
double-ring design to a single ring in order to improve the analysis of the scanning of the enriched 
BEGes. Also its position with respect to the surface was clearly stated in the specifications. 

3.2 Tests of the spectroscopic performance 

For Ovj3j3 experiments the energy resolution and the stability of the energy scale are of paramount 
importance for the identification of a peak that is expected to be rather weak. This peak, if it exists, 
is more accessible if the background can be reduced, e.g. by PSD. These arguments have motivated 
the exploration of the BEGe type detectors to substitute the standard semi-coaxial detectors. 

The second part of acceptance tests dealed with the characterization of the detector perfor- 
mance. This includes the following measurements of: (i) energy resolution of the 1332.5 keV 
gamma line of ^"Co (see Fig. ||), (ii) electronic noise component of the energy resolution (related 
to the detector capacitance) using a pulser fed to the preamplifier test input, (iii) PSD performance; 
determined with a ^^^Th gamma-ray source, using the method described in Ref. OBD and updated in 
Ref. [0] (Fig. H), (iv) long term stability of the charge collection; determined by periodic measure- 
ments with 7 sources and a pulser over an interval of several months. 

The results of this part are summarized in Tables and ^, while the measurements are de- 
scribed in the following subsections. 

3.2.1 Energy resolution 

The determination of the energy resolution follows the standard prescriptions of spectroscopy, i.e. 
to accumulate high statistic spectra with the pulser and the ^°Co source. The FWHM is the key 
parameter for the comparison. Fig. shows the zoom onto the peak at 1332.5 keV. This range is 
fitted with a step-like background and a Gaussian resulting in the widths as shown in the plots. 

The size of the steps describing the background contributions varies by a factor 2 in these five 
examples. However, due to different geometries of the detectors an interpretation of this factor is 
not meaningful, even more not important as in general this background is small, namely on the 
%-level. 

Therefore, the widths of the Gaussian can be extracted rather reliably. For comparison, the five 
Gaussians are normalized to the same height and are shown in the sixth frame on the right bottom 
of Fig. 0. The uncertainties due to the fit as shown in the plots are only a minor fraction of the total 
uncertainty. The latter is dominated by the stability of the energy scale. Therefore, the numbers 
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Figure 3. The responses of the five detectors to the 1332.5 keV 7 line of ^"Co. The spectra are fitted with a 
step-like background and a Gaussian with the FWHM as shown in the plots. The bottom right graph shows 
all Gaussians normalized to the same height for comparison. 



given in Table are reduced in the number of significant digits. All detectors except the 1/D exhibit 
a compai^able FWHM for the pulser. Similarly, 1/D has also the worst FWHM for the ^°Co line. 

The worse performance of detector 1/D is not understood. Increased pulser FWHM is an indi- 
cation of higher electronic noise, which could be e.g. caused by an increased detector capacitance, 
an effect that possibly is connected with the reduced active volume. However, a possible technical 
issue with the signal read-out electronics could not be ruled out as a source of this noise. A prob- 
lem with the charge collection as implicated in section [3.1. 6| would in general manifest itself by an 
appearance of a tail on the side of the 7 peak; however, this is not observed. 



Table 3. The energy resolution of the prototype detectors in comparison to the reference detector and the 
duration of the long-term stability tests. 



detector 


pulser 


60Co 


long-term test 


crystal/slice 


FWHM [keV] 


FWHM [keV] 


duration [d] 


reference 


0.49 


1.63 


51 


1/D 


0.74 


1.76 




2/B 


0.45 


1.60 


32 


3/B 


0.45 


1.74 


40 


4/C 


0.43 


1.62 


85 


4/D 


0.44 


1.66 


45 
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Table 4. PSD performance of the prototype detectors. To facilitate the comparison among the several 
detectors, the cut was set to reproduce a 90 % survival probability for the events in the DEP. The uncertainties 
include statistical as well as systematic contributions. The energies of the DEP, SEP, and FEP lines are given 
in keV together with the Co summation peak (SP). 



detector 


surviving fraction after PSD analysis [%] 


[crystal/slice] 


DEP 


FEPl 


SEP 


FEP2 


^^Co SP 




1592.3 


1620.6 


2103.5 


2614.5 


2505.7 


reference 


90 


10.4±0.9 


5.3±0.5 


7.7 ±0.4 


0.28 ±0.02 


1/D 


90 


14.3 ±0.5 


8.7±0.5 


13.1 ±0.5 




2/B 


90 


10.3 ±1.4 


5.0±0.6 


7.4 ±0.6 


0.093 ±0.011 


3/B 


90 


7.6±0.7 


4.3 ±0.4 


5.7 ±0.5 


0.090 ±0.012 


4/C 


90 


9.0±0.3 


4.1±0.3 


6.1 ±0.5 


0.068 ±0.011 


4/D 


90 


8.5±0.6 


4.7±0.5 


6.4 ±0.5 


0.088 ±0.009 



3.2.2 Pulse shape discrimination 

The two j3s from the Ovj3j3 event have a rather short path length in the germanium detector de- 
positing their energy in a rather small volume of a few mm^ ('single-site'). In contrast, photons of 
equivalent energy will undergo several scatterings requiring a significantly larger volume ('multi- 
site'). This behavior will be reflected in the time structure of the charge signal. To investigate the 
effect and to define analysis cuts one has to retire to 7 sources of higher energy, e.g. to ^^^Th with 
its 2614.5 keV quanta. The double escape peak (DEP) is assumed to exhibit the same behavior as 
the OVj8j8 event. The kinetic energy of the electron-positron pair is deposited in a similarly small 
volume while the two annihilation quanta escape the active volume of the detector. In contrast, 
for the full energy peak (FEP) typically involves events with multiple spatially-separated interac- 
tions, either by Compton scatterings followed by photo-absorption, or at high energies also via pair 
production and subsequent absorption of both annihilation photons. For detailed explanations and 
their validation see Refs. [0, 0]. 

The ^^^Th source is rather convenient as it provides a 1592.5 keV line, which is the DEP of 
the 2614.5 keV line in ^^^Tl In its vicinity, the FEP of the 1620.6 keV 7 of ^^^Bi appears. Both 
isotopes are ^^^Th progenies. Furthermore, both are relatively close to the Qj^p value of 2039 keV. 

The ^^^Th spectra taken with the five detectors are shown in Fig. § for the region of the DEP 
and FEP as discussed above. The original spectra are shown as red histograms. Applying the 
PSD cut after optimization and normalizing to a survival fraction of 90 % for the DEP the blue 
histograms are obtained. In order to appreciate the overall effect of PSD a large range from 600 to 
3000 keV is plotted in the frame at the bottom right of Fig. 0. 

The effect of PSD is seen clearly in the reduction of the FEP line at 1620.6 keV for all detec- 
tors. The relevant fractions are compiled in Tab. ^ where FEPl corresponds to the discussed line 
at 1620.6 keV. Except for the 1/D detector the reduction of the 'multi-site' events within the peaks 
to 8-10 % is reachable. The uncertainties include statistical as well as systematic contributions, 
the latter are derived from previous studies as described in Ref. [0, 0]. The performance is quite 
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Figure 4. Spectra of ^^^Th taken with the five detectors tested, zoomed into the range of the DEP at 
1592.3 keV and the FEP at 1620.6 keV (see text). The standard spectra are given in red (upper curves); 
applying the PSD cut reduces the data as shown in blue (lower curves). The overall effect of the PSD cut is 
demonstrated by example of detector 3/B for the range from 650 to 3000 keV (bottom, right, log scale). 



similar for all other detectors. 

The bottom right panel in Fig. ^ showing the complete spectrum above 650 keV demonstrates 
the effect of the PSD cut for a wide energy range. Actually, it hints at the fraction of 'multi-site' 
events contributing to the total spectrum. The separation of SSE and MSB is energy independent 
which was validated in the range from 1.4 to 2.4 MeV, interesting for OVj8j3 decay. The reduction 
is more or less constant at 30 % up to the Compton edge around 2300 keV, only above that energy 
the intensity is reduced to about 5 % by the PSD cut due to the prevalence of MSB. 

Additional information is provided for the single escape peak (SBP) at 2103.5 keV and the 
2614.5 keV line of ^^^Tl (FBP2). The last column of Tab. | shows the PSD performance for the 
summation peak (SP) of the two ^"Co lines that is seen at an energy of 2505.7 keV. The summa- 
tion effect results in highly 'multi-site' events, responsible for more than two orders of magnitude 
stronger suppression than is reached in the FBP and SBP. The suppression achieved with the four 
tested detectors is in this case much better than that quoted for the reference BBGe, as the analysis 
procedure has ben improved meanwhile. 

The PSD performance of the detector 1/D was again anomalous with higher than typical back- 
ground survival fractions. The cause could possibly be due to the higher electronic noise present 
in the measurements of this detector, or due to less favorable electric field distribution inside the 
detector, possibly related to the charge collection deficiency indicated in section [3.1.6| . 
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Figure 5. Long term stability test of detector 4/C with two sources. The blue crosses represent the ^"Co 
data whereas the red diamonds show the results for Am. 



3.2.3 Long-term stability 

The last column of Table gives the duration of the long term tests lasting between 1 and 3 months. 
The diodes were iiTadiated alternating with an Am and a ^"Co source. The count rates were nor- 
malized to the start rate and plotted in Fig. ^ together with the stability of peak position and energy 
FWHM. The count rate stays well within ±1 %. Measurement of possible charge collection effects 
on peak position variation was limited by the inherent variability of the DAQ system gain (charac- 
terized via shorter pulser stabiUty measurements) to the upper limit of about ±0.01 % (±0.1 keV) 
at the 1332.5 keV line and about ±0.05 % (±0.05 keV) at the 59.5 keV line. The FWHM re- 
mained also very stable within approximately ±20 eV. Similar results were obtained with the other 
detectors, except of 1/D, for which a long term test was not performed. 

4. Conclusions 

The whole production chain, from material procurement to detector operation, was successfully 
tested with depleted germanium left over after the enrichment process. Various detectors were 
produced with modifications of the production processes with the purpose to maximize the yield 
of the upcoming new BEGe detectors from the available enriched Ge. The latter, after reduction 
and purification at PPM Pure Metals, amounts to 35.4 kg plus a 1.1 kg tail of slightly lower than 
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6N quality which can be recycled in a subsequent production run. The production of new BEGe 
detectors for the second phase of Gerda and their characterization is ongoing. 

The measured characteristics and performance of the tested depleted detectors were as good 
as those of a reference BEGe from standard Canberra production made of natural germanium. This 
proves that high-purity germanium material crystals produced from isotopic modified germanium 
raw material according to the described production steps have the same properties as those made 
from standard material and processing by Canberra. The depleted germanium has the same chemi- 
cal history as the enriched germanium for the second phase of Gerda. Thus, this test confirms the 
viability of the production chain for enriched BEGe detector production. 

Good-performance detectors were produced successfully from germanium crystal slices which 
are normally discarded before BEGe production. In addition, the surfaces of some slices were 
machined less than usual with no loss of performance. Thus, the smaller material losses will 
increase the yield of total detector mass from the limited available enriched germanium. Presently, 
the yield of detector mass was about 50 % with respect to the 6N starting material. It can be 
further augmented by increasing the initially available enriched germanium, as in the case of the 
subsequent enriched germanium production. The residual material, including seed and tail ends of 
crystals, as well as kerf and grinding losses can be recovered and chemically reprocessed, and a 
large fraction can return into the crystal production process. 

The viability of the production processes of isotopically modified BEGe detectors for the 
OVj8j8 decay experiment Gerda was proven within a rather short period of less than 2 years in- 
cluding all steps. The vital part for the enriched program amounted to about 60 % of this time. 

A precise positioning system for the radioactive sources with well known intensity is of 
paramount importance for the characterization process of new detectors. Such a system was es- 
tablished for the tests of the enriched detectors. For a Ovj3j3 experiment the precision of the active 
volume is the parameter of highest relevance followed by the energy resolution itself. 
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